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Prediction of the derivative discontinuity in density functional theory from an 
electrostatic description of the exchange and correlation potential 

Xavier Andrade^'jj and Alan Aspuru-Guzik^'jj 

^Department of Chemistry and Chemical Biology, Harvard University, 
12 Oxford Street, Cambridge, MA 02138, United States 

We propose a new approach to approximate the exchange and correlation (XC) functional in 
density functional theory. The XC potential is considered as an electrostatic potential, generated 
by a fictitious XC density, which is in turn a functional of the electronic density. We apply the 
approach to develop a correction scheme that fixes the asymptotic behavior of any approximated XC 
potential for finite systems. Additionally, the correction procedure gives the value of the derivative 
discontinuity; therefore it can directly predict the fundamental gap as a ground-state property. 
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An important and long standing topic in density func- 
tional theory (DFT) U is the prediction of the funda- 
mental gap [2 [3] -Eg, which is defined as the difference of 
the ionization energy and the electron affinity. In DFT, 
the gap is not simply the difference between the Kohn- 
Sham (KS) eigenvalues of the highest occupied molecular 
orbital (HOMO) and the lowest unoccupied molecular or- 
bital (LUMO). Instead, it is given by ||] 



E^ — clumo — chomo + Ax 



(1) 



where chomo and eLUMO are the HOMO and LUMO 
KS eigenvalues, respectively, and Axe is the derivative 
discontinuity (DD) of the XC energy with respect to the 
particle number N, 



Axe = 14e(iV+) - V^,iN') 



(2) 



For the local density approximation (LDA) and many 
generalized gradient approximations (GGA), the DD is 
zero [SI . In these approximations the predicted gap is ef- 
fectively the KS gap, which severely underestimates the 
experimental value. Even for functionals that are discon- 
tinuous with the particle number, the DD is not simple 
to calculate [31 |BHH] ■ Alternatively, XC approximations 
have been proposed where the KS gap is directly used to 
predict the gap [9^ avoiding the calculation of the DD. 

In this article, we present an XC functional for finite 
systems that, with similar computational cost as a LDA 
or GGA calculation, has the right asymptotic limit for 
low density regions and directly provides the value of 
the DD. Hence, the proposed functional can predict the 
fundamental gap for an atom or molecule as a ground 
state property. 

The approach that we advocate is not based on increas- 
ing the number of functional variables, but on changing 
the way that the XC potential is described: we consider 
the XC potential as an electrostatic potential, generated 
by a fictitious XC charge density. In contrast to di- 
rectly modelling the potential, the XC density becomes 
the quantity to approximate as a functional of the elec- 
tronic density n. 



Given a XC potential V^c we define the XC density n^xc 
by the Poisson equation (atomic units are used through- 
out) 



Wxe(r) = -47rnxe(T-) 



(3) 



with the boundary condition V^df — ^ oo) = 0. 

The justification for this approach comes from an 
important property of the XC potential: the so-called 
asymptotic limit [TUl E] , 



Ke(r) 



\r\ ^oo) , 



(4) 



The LDA and most GGAs do not obey the asymptotic 
limit condition. In part, this common deficiency can be 
explained simply. In regions that are spatially far away 
from the system, the density and its derivatives decay 
exponentially to zero [H] . It is difficult to use local values 
of the density and its derivatives to reproduce a field that 
decays to zero much more slowly. 

For the XC density the asymptotic limit implies two 
simple conditions: normalization and localization. For- 
mally, 

n^c{r) = (|r| -^ cx)) , (5a) 

/"drnxe(r) = -l , (5b) 

(proof in supp. mat. |12j). These constraints are simi- 
lar to the ones for the electronic density. So, in princi- 
ple, it is simple to construct a local or semi-local density 
functional nxc[?T-] with the proper asymptotic limit. In 
fact, a direct example of this XC density based approach 
is the functional n^c[n\{r) = —n{r)/N, which yields the 
Fermi- Amaldi XC potential [13] and provides an accurate 
approximation of V^c in the asymptotic regime [141 115] . 
It is illustrative to see what the XC density looks like 
for standard DFT functionals. For a given XC potential, 
rixe can be calculated analytically (see supp. mat. 12 ), 
but in practice it is simpler to evaluate Eq. ([3| numer- 
ically. In Fig. [T] we compare nxc for LDA and exact 
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FIG. 1. KLI electronic density and KLI and LDA XC densi- 
ties for the beryllium atom. The inset details the asymptotic 
region, where the KLI XC density shows the correct behav- 
ior. The LDA XC density is incorrect in this region, with a 
positive part that screens the charge in the central region. 



exchange in the Krieger-Li-Iafrate (KLI) [15] approxima- 
tion. As can be seen from the figure, for large radius the 
KLI XC density correctly goes to zero faster than the 
electronic density. On the other hand, the LDA XC den- 
sity becomes positive. This positive XC charge screens 
the XC charge in the central region making the total XC 
charge zero and therefore causing the potential to decay 
exponentially. 

As a first application of our approach, we propose a 
correction method to enforce the proper asymptotic limit 
for functionals that do not have it by construction. Given 
a certain potential V^c, we calculate the associated n-^c 
from Eq. ((sl) . To this XC density we apply the correction 
procedure, that generates a corrected XC density n^^,. In 
turn, the corrected XC density is used to reconstruct the 
corrected XC potential V^^ by solving Eq. (|3|. 

The correction procedure for fixe enforces it to be local- 
ized by setting it to zero when the local value of electronic 
density is below a certain threshold rj. This simple pro- 
cedure can be written as a correction term Anxc to be 
added to the XC density of the original functional. 



AnxcM(r,77) = 



if n{r) > rj 

-nxcN(r) if n(r) < 77 



(6) 



To determine the parameter 77, for each density we ob- 



tain an optimized value 770 that tries to enforce Eq. ( 5b I 



First, we define the total XC charge as a function of ?/ 



&c(»7) 



/ dr {rix 



M(r) + AnxcM(r,ry)} . (7) 



Ideally, from Eq. (5b), we need to find 770 such that 
Qxcivo) — ^1- However, there is no guarantee about the 
existence or uniqueness of rjQ. So we choose t/q such that 
Qxcivo) has the closest value to -1, with 7^0 restricted to 
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FIG. 2. CXD-LDA potential for neon. Comparison with the 
LDA, LB [m, R.PP [20] and the exact [M] potentials. 



be smaller than the first minimum of (7xc('?)- (See supp. 
mat. [H].) 
When (7x( 



^(t^o) 7^ —1, Eq. (5b) is still not satisfied, so we 
rescale the correction by |(7xc(??o)| ■ The final expression 
for the XC density of the corrected functional is 

1 



.(r) 



.(r) 



kxc(f7o)| 



A77.xc(r,7/o) 



(8) 



This rescaling form guarantees that Eq. (5b) is satisfied. 



and that the original XC potential is changed as little as 
possible in the central region (the region where n > rjo). 

In theory, it is only the exchange term that is responsi- 
ble for the long range behavior, due to the much faster de- 
cay of the correlation term [11], therefore, the correction 
can be applied either to the exchange potential or the full 
XC potential. In this work, we apply it to the exchange 
part of the LDA functional, and we call the combination 
of the corrected LDA exchange and LDA correlation (in 
the Perdew- Wang form [17] ) the corrected exchange den- 
sity LDA (CXD-LDA). For spin-polarized systems, the 
correction is calculated for the spin-unpolarized poten- 
tial using the total density. Then the difference between 
the corrected potential and the original one is added to 
the XC potential for each spin component. 

To test the CXD-LDA functional, we performed cal- 
culations for atoms, hydrogen to strontium, and for a 
set of small molecules. We implemented the correction 
procedure in the APE [TH] and Octopus [IH] codes. We 
find that the optimized value of rjo changes significantly 
for different systems. For most atoms gxc(?yo) 7^ —1 
while for all the tested molecules gxc(??o) = — 1 (see supp. 
mat. [H]). The numerical cost of a self-consistent solu- 
tion using the correction is similar to the LDA calculation 
(see supp. mat. for details [T2]). 

In Fig.|2j we show the CXD-LDA potential for Ne com- 
pared with an accurate approximation to the exact po- 
tential [Uj, the original LDA functional and two func- 
tionals that have the correct asymptotic behavior: the 



van Leeuwen-Baerends (LB) GGA [TT] and the Rasanen- 
Pittalis-Proetto (RPP) meta-GGA [20]. 

For atoms it is simple to understand the effect of the 
correction procedure. Due to the spherical symmetry 
and the monotonically-decreasing density, the correction 
XC charge Arixc is a spherical shell. By Newton's shell 
theorem, the correction potential will be constant in the 
central region. Outside, the correction will decay close 
to — 1/r. We can expect this behavior to be similar for 
more complex systems if the n(r) = 7/0 surface is close to 
a sphere (see supp. mat. [H]). 

The shape of the correction is similar to the one pro- 
posed by Casida and Salahub [211 [22], who argument that 
a shift of the XC potential in the central region is nec- 
essary to fix the asymptotic limit of the LDA potential. 
Moreover, they show that the shift is related to the DD 
of the energy with respect to the particle number. As 
in our method the shift appears naturally from imposing 
the asymptotic limit, we can obtain the value of the DD. 

To obtain the relation between the DD and the shift, 
we assume that a potential V^c, which lacks the DD, ap- 
proximates the XC potential averaged over the disconti- 
nuity [SIT]. This is 



K< 



[V^ciN+) + V^ciN-)] 



Using Eq. (pi), immediately follows that 



V^c{N-) = Kc - 2^- 



(9) 



(10) 



By imposing the asymptotic limit of Eq. |4l our cor- 
rected potential is approximating Vxc{N~) [23]. There- 



fore, we can obtain the value of the DD from Eq. ( 10 1 



For practical calculations we average the change of the 
XC potential due to the correction over the central region 



A. 



/ dr[K'c(r-)-Kc(r)] , (11) 



where fi is the volume of the central region. This ex- 
pression for the DD can be calculated directly from the 
correction process as a ground-state property. An alter- 
native, but less practical, method for the calculation of 
the DD is detailed in supp. mat. [12j. 

In Table llj we compare the DD obtained with Eq. ( 11 ) 
with the values reported by Chan fS^ from ensemble DFT 
(with XC potentials obtained from wave-function meth- 
ods). We also compare it with the experimental value 
of the gap, that for these open-shell atoms is equal to 
the DD since the KS gap is zero. The three sets present 
a remarkable agreement, with our results being smaller 
that the experimental values by less than 10%. 

To investigate further the quality of corrected func- 
tional and its DD, we compare the calculated gap with 
the LDA KS gap and the experimental gap, for our set of 
atoms and molecules. The results are plotted in Fig. [3J 



TABLE I. Comparison of the calculated derivative discontinu- 
ity with theoretical and experimental results. " Experimental 
gap values from Ref. [24]. *" Ensemble spin DFT results by 
Chan [6]. Values in atomic units. 

Atom Experimental" ESDFT'' CXD-LDA 
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0.295 


0.270 


0.284 
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0.367 


0.342 


0.337 
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FIG. 3. Comparison of the CXD-LDA gap with LDA and 
available experimental results. The CXD-LDA gap has con- 
tributions from the Kohn-Sham (KS) gap and the derivative 
discontinuity (DD). Result for atoms (top) from H to Kr, 
experimental values from Ref. (2^, and a set of molecules 
(bottom), experimental values compiled in Ref. [7], except 
for CeHe [HHl]. (Data in supp. mat. [H].) 



The KS gap of the corrected functional is close to the 
LDA one [2F and far from the experimental value. Once 
we add the DD, however, the results are closer to the 
experiment, with an average error of 11% for atoms and 
7% for molecules. 

While the correction has little effect on the KS gap, it 
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FIG. 4. LDA and CXD-LDA ionization energy as a function 
of the experimental value. Atoms from H to Sr. (Data in 
supp. mat. [E].) 



TABLE IL Mean absolute error in the ionization energy for 
atoms (He to Ar). Comparison of CXD-LDA with LDA, 
LB [TT], RPP Uni and KLLCS [ISlllZl functionals. "Results 
from ReL [28]. 'Results from Ref. [29] given in Ref. [28] . 
(Data 



in supp. mat. 


|12|.) 






LDA LB" 


RPP" 


KLLCS'' 


CXD-LDA 


41% 3.7% 


7.4% 


5.7% 


4.1% 



changes the KS eigenvalues. This can be seen in the ion- 
ization energy (I), which in DFT is given by — enoMO [ID]- 
In Fig.|4) we plot I for the LDA and CXD-LDA function- 
als as a function of the experimental value. In Table [llj 
we compare the deviation from experimental results for 
atoms with other XC functionals that have the proper 
asymptotic limit: LB, RPP, and KLI with Colle-Salvetti 
correlation [27 (KLI-CS). The correction procedure im- 
proves considerably the LDA results, with similar accu- 
racy to other long range XC potentials. 

In summary, we have introduced a new auxiliary quan- 
tity, the XC density, to construct approximations for the 
XC potential. Based on an exact condition that the XC 
potential must fulfill and basic notions of electrostatics, 
we have presented a correction method for any previously 
proposed XC potential. 

Additionally, the correction procedure allows for the 
direct calculation of the DD of the XC energy, which 
can be used to directly predict the fundamental gap as 
a ground-state property. Moreover, our approach al- 
lows for a routine computation of the DD as a practical 
method for the prediction of the gap. 

The proposed potential is a pure functional of the elec- 
tronic density with a certain degree of non-locality in- 
cluded by the optimization of r?o and by the Poisson equa- 
tion. The correction procedure does not depend on any 
empirical or globally adjusted parameter. 



Since the basis for our method is the correction of the 
XC potential in the asymptotic region, it is not directly 
applicable to crystalline systems. However, the concepts 
of the XC density and the DD as a potential shift are still 
valid. Therefore, it might be possible to generalize the 
method to solids, where the determination of accurate 
gaps is one of the main challenges for DFT. 

We acknowledge support from the US National Science 
Foundation (DMR-0934480). 
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Supplementary Material 
Analytic expression for the XC density 

Given an XC potential, the XC density is given by 

1 



by using the chain rule for functional derivatives 



(12) 



n^c[n]ir) 



1 

An 



Jdr'U,[n](r,r')\/Mr') 
/"dr'dr"fcxcW(r,r',r") Vn(r') • Vn(r") 



(13) 



where /xc and fcxc are, respectively, the first and sec- 
ond functional derivatives of V^c with respect to the den- 



sity. Interestingly, Eqs. (12) and (13) give a recipe to 



reconstruct an XC potential from its first two functional 
derivatives. 



Proof of the conditions for the XC density 

The asymptotic condition for V^c implies that there 
exists a certain r^ such that 



Kc(r) 



1 



\r\ > r,) . 



(14) 



If we apply the Laplacian over this equality, we get 

nxc(r) = (|r| > r,) , (15) 

the localization condition. 



For the normalization condition, we integrate Eq. ( 12 ) 



over a spherical volume of radius r^ and we use Gauss's 
theorem over the right-hand side 



drnxc(r) = --— 
J, >r 47r 



dA-VKcc(r) 



(16) 



From Eq. ( 14 ), V14c["-] = r/r'^ which is constant over the 



sphere, so the integral on the right is just l/r^ times the 
surface of the sphere, Airr"^. The integral on the left hand 



size can be extended to whole space due to Eq. (15), so 
finally 



/ 



drn^^ir) 



-1 . 



(17) 



Optimization of the rj parameter 

Given a XC density rixc, the corrected XC density is 
defined as 



7^xc(T■) = nxc(r) + 



1 



lfec(?7)l 



Anxc(r-,r?), (18) 




FIG. 5. The q^dv) function, Eq. ( |20| ), for LDA on atomic 
Beryllium atom and CeHe. The crosses mark the position of 
•qo, the optimal value of 77. 



with the correction XC density 

[ -~n^c[n\[r) li n{r) < 77 
and the total XC charge 

q^civ) = / dr {nxcW(r) -I- AnxcM(r,r/)} . (20) 

The value of rj is optimized for each density, the optimum 
value 770 is the value that has the total charge closest to 
-1. Formally, 



r]o[n] = min |gxcN(?/) 



II 



(21) 



where r]rn is the smallest value where (7xc('7) has a mini- 
mum. These conditions are necessary, since it is not al- 
ways possible to find a unique value where qxdv) = ^1- 
This can be seen in Fig. ^were we show q^dv) for two 
systems. While for beryllium ^xc never reaches -1, for 
benzene there are two points where it does. The marks 
in the figure indicate the value of ryo and (7xc('7o) selected 
by our method. 



Corrected functional test 

To test the corrected LDA functional, CXD-LDA, we 
did calculations for atomic systems from hydrogen to 
strontium and for a set of molecules. 

The optimized parameter 770 obtained for each system 
is plotted Figs. [6] and Figs.[7| For atoms we also plot the 
corresponding value of qxciva) and Vc, the value where 
n{rc) = rjo, that marks the position where the XC den- 
sity is set to zero. For all tested molecules q^dvo) = ~1- 
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FIG. 6. CXD-LDA optimized parameters for atoms. Top: the 
density threshold rio- Center: the radius re where n{rc) = rjo. 
Bottom: The value of the optimized total XC charge, (j'xc(?7o)- 
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FIG. 8. CXD-LDA potential for helium. Comparison with 
the LDA, LB [11], RPP [SO] and the exact [M] potentials. 
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FIG. 7. Optimized density threshold, rjo, for a set of small 
molecules. The value of the total XC charge, gxc(??o), is -1 for 
all cases. Values are given in Table |IV| 



In Fig. |8j we plot the XC potential for He compared with 
other functionals. In Fig. 10 we present a more detailed 
version of the plot for the band gap of the atoms. Finally, 
we give the values for our optimized parameters, ioniza- 
tion energies and gaps in Tables IV (molecules) and [V| 
(atoms). 



Shape of the correction 

In Fig. [9] we plot the difference of the corrected XC po- 
tential and the original one for 5 atoms and 2 molecules. 
For atoms the correction has the shape expected from 
electrostatic principles. For molecules, the correction in 
the central region is not exactly constant. This is ex- 
pected as the n — tjq surface is not necessarily a perfect 
sphere. However the deviation from a constant shift is 
not large and we can still obtain a representative value 



from the average. Near the atomic positions there are 
some peaks due to numerical error in the finite difference 
calculation of the Laplacian, but these are a few points 
that do not influence considerably the average. 



Alternative calculation of the derivative 
discontinuity 

We propose to calculate the derivative discontinuity 
(DD) from the shift of the potential due to the correction 
procedure 



A^ =- 



n 



n(r)>?7 



dr[l4c(r)-yxc(r)] 



(22) 



Alternatively, the DD can be obtained by comparing the 
results of the corrected XC potential with the results of 
the uncorrected functional. In principle, in DFT the ion- 
ization energy is given by the highest occupied Kohn- 
Sham eigenvalue 



ehomo 



(23) 



however, for a potential that averages over the disconti- 
nuity [7] 



ehomo 



-/ 



1 



:A. 



So, the DD is also given by 

^xc — 2(eHOMO ^ chomo) • 



(24) 



(25) 



In table |III| we give the values of the DD obtained 



with the two previously mentioned approaches, Eq. (22) 



and Eqs. (23) and (24), for four different atoms. 
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FIG. 9. Correction of the LDA X potential for different atoms 
(top) and molecules (bottom). The vertical lines mark the 
position where n = rjQ. For the HCN molecule the value of 
the shift is plotted for the axis that passes by the atoms (HCN 
par.) and other axis that is perpendicular to that one (HCN 
perp.). For CeHe one axis in the plane of the molecule (CeHg 
in plane) and the other is perpendicular to that plane (CeHg 
perp.). 



B 


0.281 


0.284 


C 


0.333 


0.337 





0.427 


0.435 


F 


0.460 


0.467 



TABLE HI. Comparison of methods to calculate the deriva- 
tive discontinuity, Eq. (22 1 and Eq. (25 I. 



Computational cost of the correction 



In general the cost of the correction procedure is not 
significant in comparison to the total cost of a DFT self- 
consistent procedure. The additional cost conies from 



the numerical solution of Eq. ( 12 ). This is the same pro- 



cedure required to obtain the Hartree potential and can 
be done efficiently in linear or quasi-linear computational 
time using methods like fast Fourier transforms, multi- 
grid or fast multipole expansions. Moreover, due to lin- 
earity, it would be possible to solve the Poisson equation 
once, obtaining at the same time the Hartree and the 
corrected XC potentials. 
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FIG. 10. Comparison of fundamental gap calculations with available LDA and CXD-LDA functionals with experimental results 
for atoms between H and Sr. Experimental values obtained from Ref. [23]. Values are given in Table Iv] 



Molecule 


Experimentar 

E, 


L 
Aks 
0.236 


DA 

AEg % 

66 


no 

1.5E-03 


Ak. 
0.249 


CXD-LDA 

Axe 
0.356 


E, 
0.605 


AEg % 


H2O 


0.70 


13.6 


NH3 


0.60 


0.204 


66 


6.4E-04 


0.214 


0.317 


0.532 


11.4 


CH4 


0.79 


0.335 


58 


2.5E-04 


0.348 


0.294 


0.642 


18.8 


HP 


0.81 


0.322 


60 


4.1E-03 


0.341 


0.413 


0.754 


6.9 


PH3 


0.46 


0.220 


52 


1.3E-04 


0.225 


0.256 


0.482 


4.7 


CO 


0.58 


0.245 


58 


3.8E-04 


0.244 


0.312 


0.556 


4.1 


Na 


0.65 


0.286 


56 


3.7E-04 


0.286 


0.318 


0.604 


7.1 


HaS 


0.46 


0.210 


54 


2.3E-04 


0.210 


0.274 


0.484 


5.2 


HCN 


0.58 


0.281 


52 


2.4E-04 


0.282 


0.292 


0.574 


1.1 


C2H2 


0.52 


0.245 


53 


1.7E-04 


0.245 


0.270 


0.515 


0.9 


C2H4 


0.46 


0.207 


55 


8.9E-05 


0.208 


0.255 


0.463 


0.6 


H2CO 


0.46 


0.122 


74 


1.6E-04 


0.121 


0.287 


0.408 


11.3 


CgHg 


0.38 


0.188 


51 


1.3E-05 


0.188 


0.199 


0.386 


1.4 



TABLE IV. Experimental, LDA and CXD-LDA results for a set of molecules, ryo: optimized density threshold value from 
Eq. (|19[), in all cases qxc(7?o) = — 1. Aks: HOMO-LUMO gap, Axe: derivative discontinuity. Eg-, fundamental gap, for CXD- 



LDA is Aks + Axe, and AEg-. Difference in the gap with experimental results, for LDA Aks is considered as the gap. All 
calculation are done with the Octopus code 19 using a grid composed of spheres of radius 16 a.u. around each atom with 
uniform spacing 0.25 a.u, and LDA norm-conserving pseudo-potentials. Geometries from Ref. [30^. " Experimental values 
compiled in Ref. [7j, except for CeHe obtained as difference between the ionization energy and the electron affinity from 
Refs. [MESl. 
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Atom 


Experimentar 

I Eg 


^horno 

-0.269 


L 
A/% 
46 


DA 

Aks 
0.173 


AEg % 

63 


Vo 


9xc(??o) 


^hoino 


CXD-I 

A7% 


.DA 

Aks 


Axe 


E, 


AEg % 


H 


0.500 


0.472 


6.8E-04 


-0.53 


-0.421 


15.8 


0.174 


0.307 


0.481 


2.0 


He 


0.904 


- 


-0.570 


37 


0.570 


- 


8.2E-03 


-0.63 


-0.804 


11.0 


0.655 


0.470 


1.125 


- 


Li 


0.198 


0.175 


-0.116 


41 


0.042 


76 


l.lE-04 


-0.69 


-0.205 


3.5 


0.043 


0.177 


0.220 


25.4 


Be 


0.343 


- 


-0.206 


40 


0.129 


- 


7.1E-04 


-0.82 


-0.326 


4.8 


0.129 


0.241 


0.371 


- 


B 


0.305 


0.295 


-0.151 


51 


0.000 


100 


2.5E-03 


-0.86 


-0.291 


4.5 


0.000 


0.284 


0.284 


3.5 


C 


0.414 


0.367 


-0.227 


45 


0.000 


100 


6.2E-03 


-0.92 


-0.394 


4.8 


0.000 


0.337 


0.337 


8.2 


N 


0.534 


0.534 


-0.308 


42 


0.148 


72 


l.lE-02 


-0.96 


-0.502 


6.1 


0.148 


0.391 


0.539 


0.9 


O 


0.501 


0.447 


-0.272 


46 


0.000 


100 


2.2E-02 


-1.00 


-0.486 


2.9 


0.000 


0.435 


0.435 


2.7 


F 


0.640 


0.515 


-0.384 


40 


0.000 


100 


l.lE-02 


-1.00 


-0.614 


4.1 


0.000 


0.467 


0.467 


9.4 


Ne 


0.792 


- 


-0.498 


37 


0.498 


- 


l.OE-02 


-1.00 


-0.746 


5.9 


0.560 


0.501 


1.061 


- 


Na 


0.189 


0.169 


-0.113 


40 


0.034 


80 


l.lE-04 


-0.69 


-0.198 


4.7 


0.034 


0.169 


0.203 


20.5 


Mg 


0.281 


- 


-0.175 


38 


0.125 


- 


4.7E-04 


-0.86 


-0.282 


0.5 


0.126 


0.215 


0.341 


- 


Al 


0.220 


0.204 


-0.111 


50 


0.000 


100 


1.3E-03 


-0.91 


-0.222 


0.8 


0.000 


0.224 


0.224 


9.8 


Si 


0.300 


0.249 


-0.170 


43 


0.000 


100 


1.2E-03 


-1.00 


-0.297 


1.0 


0.000 


0.256 


0.256 


2.9 


P 


0.385 


0.358 


-0.231 


40 


0.079 


78 


6.8E-04 


-1.00 


-0.368 


4.5 


0.080 


0.275 


0.355 


0.9 


S 


0.381 


0.304 


-0.229 


40 


0.000 


100 


7.3E-04 


-1.00 


-0.375 


1.6 


0.000 


0.293 


0.293 


3.6 


CI 


0.477 


0.344 


-0.305 


36 


0.000 


100 


7.3E-04 


-1.00 


-0.462 


3.0 


0.000 


0.315 


0.315 


8.5 


Ar 


0.579 


- 


-0.382 


34 


0.373 


- 


7.4E-04 


-1.00 


-0.549 


5.2 


0.394 


0.335 


0.729 


- 


K 


0.160 


0.141 


-0.096 


40 


0.023 


84 


7.5E-05 


-0.74 


-0.171 


6.9 


0.023 


0.149 


0.172 


22.2 


Ca 


0.225 


0.224 


-0.142 


37 


0.088 


61 


2.7E-04 


-0.94 


-0.232 


3.4 


0.089 


0.182 


0.271 


21.0 


Sc 


0.241 


0.234 


-0.141 


42 


0.000 


100 


4.0E-04 


-0.95 


-0.238 


1.1 


0.000 


0.195 


0.195 


16.6 


Ti 


0.251 


0.248 


-0.175 


30 


0.000 


100 


5.3E-04 


-0.95 


-0.277 


10.2 


0.000 


0.205 


0.205 


17.3 


V 


0.248 


0.229 


-0.186 


25 


0.000 


100 


6.4E-04 


-0.94 


-0.292 


17.8 


0.000 


0.212 


0.212 


7.1 


Cr 


0.249 


0.224 


-0.196 


21 


0.000 


100 


7.0E-04 


-0.93 


-0.305 


22.7 


0.000 


0.219 


0.219 


2.3 


Mn 


0.273 


- 


-0.206 


25 


0.045 


- 


8.4E-04 


-0.92 


-0.317 


16.1 


0.045 


0.224 


0.269 


- 


Fe 


0.290 


0.285 


-0.183 


37 


0.000 


100 


9.5E-04 


-0.92 


-0.299 


2.9 


0.000 


0.233 


0.233 


18.2 


Co 


0.290 


0.265 


-0.194 


33 


0.000 


100 


1.2E-03 


-0.92 


-0.314 


8.5 


0.000 


0.241 


0.241 


9.2 


Ni 


0.281 


0.238 


-0.205 


27 


0.000 


100 


1.3E-03 


-0.92 


-0.328 


16.7 


0.000 


0.248 


0.248 


4.0 


Cu 


0.284 


0.239 


-0.184 


35 


0.030 


88 


3.4E-03 


-0.83 


-0.309 


8.8 


0.030 


0.253 


0.283 


18.5 


Zn 


0.345 


- 


-0.223 


35 


0.176 


- 


1.7E-03 


-0.91 


-0.351 


1.8 


0.178 


0.259 


0.437 


- 


Ga 


0.220 


0.205 


-0.110 


50 


0.000 


100 


2.5E-03 


-0.91 


-0.226 


2.3 


0.000 


0.236 


0.236 


15.2 


Ge 


0.290 


0.245 


-0.165 


43 


0.000 


100 


1.6E-03 


-1.00 


-0.291 


0.1 


0.000 


0.253 


0.253 


3.2 


As 


0.361 


0.331 


-0.220 


39 


0.070 


79 


4.9E-04 


-1.00 


-0.351 


2.7 


0.070 


0.262 


0.332 


0.5 


Se 


0.358 


0.284 


-0.218 


39 


0.000 




4.6E-04 


-1.00 


-0.355 


1.0 


0.000 


0.272 


0.272 


4.3 


Br 


0.434 


0.311 


-0.283 


35 


0.000 


100 


3.9E-04 


-0.99 


-0.425 


2.1 


0.000 


0.285 


0.285 


8.2 


Kr 


0.514 


- 


-0.346 


33 


0.332 


- 


3.6E-04 


-1.00 


-0.497 


3.4 


0.347 


0.302 


0.649 


- 


Rb 


0.154 


0.136 


-0.092 


40 


0.020 


85 


6.6E-05 


-0.76 


-0.164 


6.7 


0.020 


0.144 


0.164 


21.0 


Sr 


0.209 


0.207 


-0.132 


37 


0.070 


66 


2.5E-04 


-0.96 


-0.218 


4.0 


0.070 


0.172 


0.242 


16.9 



TABLE V. Experimental, LDA and CXD-LDA results for atomic systems. 770: optimized density threshold value from Eq. ( 19 1, 
qxc('7o): optimized value of the total XC charge [Eq. (20l], ehomo: HOMO eigenvalue, /: Ionization energy, A7: difference 
between — ehomo and the experimental /, Aks: HOMO-LUMO gap. Axe: derivative discontinuity. Eg-, fundamental gap, for 
CXD-LDA is Aks -|- Axe, and AEg-. difference in the gap with experimental results, for LDA Aks is considered as the gap. 
Values calculated using the APE code [18] with spin-polarization. " Experimental results from Ref. [24], Eg values calculated 
as the difference of the ionization energy and the electron affinity. 



